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^ ■ Abstract 

D ■ Renormdynamic equations of motion and their solutions are given. New equation 

Q . for NBD distribution and Riemann zeta function invented. Explicit forms of the 

'sf I z-Scaling functions are constructed. 
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^ 1 Renormdynamics 



C/3 



In quantum field theory (QFT) existence of a given theory means, that we can control 
its behavior at some scales (short or large distances) by renormalization theory [H [2] ■ If the 
theory exists, than we want to solve it, which means to determine what happens on other 
^1 (large or short) scales. This is the problem (and content) of Renormdynamics. The result 
r-| I of the Renormdynamics, the solution of its discrete or continual motion equations, is the 
effective QFT on a given scale (different from the initial one). We can invent scale variable A 
and consider QFT on D + 1 + 1 dimensional space-time-scale. For the scale variable A G (0, 1] 
it is natural to consider g- discretization, < g < 1, A„ = g", n = 0, 1,2, ... and p - adic, 
Q> I nonarchimedian metric, with = P - prime integer number. The field variable (p{x,t,X) 
^ I is complex function of the real, x, t, and p - adic. A, variables. The solution of the UV 
1/-^ ■ renormdynamic problem means, to find evolution from finite to small scales with respect to 
the scale (time) r = InA/Ao G (0,— oo). Solution of the IR renormdynamic problem means 
to find evolution from finite to the large scales, r = InA/Ao G (0, oo). This evolution is 
determined by Renormdynamic motion equations with respect to the scale-time. 
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'^I 1.1 Renormdynamics of Q CD 
^ ! 

i The Renormdynamic (RD) equations play an important role in our understanding of Quan- 
tum Chromodynamics (QCD) and the strong interactions. The beta function is the most 
prominent object for QCD RD equations. The calculation of the one-loop /3-function in QCD 
has lead to the discovery of asymptotic freedom in this model and to the establishment of 
QCD as the theory of strong interactions [21 HI [5j. 

The MS-scheme [6] belongs to the class of massless schemes where the /3-function does not 
depend on masses of the theory and the first two coefficients of the /3-function are scheme- 
independent. 

RD equation for the coupling constant of QCD belongs to the following class of equations, 

d = /3ia + /32a^ + ... = (1) 



n>l 



where, in the case of the QCD, 



; A)2 ^ = ^LA^ = - —Nc, ... (2) 
'An' 2 ' 3 3 ^ ' 



The equation ([T]) can be reparametrized, 

ait) = fiAit)) =A + hA' + ... + UA" + ... = E /-^^ 

n>l 

A = b,A + b^A^ + ... = E (3) 
a = Af'iA) = {b,A + b^A' + ...)(1 + + ... + + ...) 

= MA + f2A^ + ... + fnA" + ...) + (32{A^ + 2hA' + ...) + ... + /3„(A" + nhA'^^^ + ...) + ... 
= + (/32 + /3i/2)A' + (/33 + 2/32/2 + Pih)A' + ... + + (n - l)/3„_i/2 + ... + + ... 

= E ^"&nin2/n25n,ni+n2-l = E PrnfT^ ■ ■ J^'' f {n, m.mi, .. . ,mk) , 

n,ni,?i2>l n,m>l;mi ,...,mfc>0 

/(n, m, mi, mfc) = ; r5n,mi+2m2 + ...+fcmi,5m,mi+m2 + ...+mfc) (4) 

mi!...mfc! 

&1 = A, &2 = /32 + /2/3l - 2/2&I = /32 - /2/3l, 

63 = /33 + 2/2/32 + Ml - 2/2&2 - 3/361 = /33 + 2(/2 - /3)/3i, 

64 = /34 + 3/2/33 + /22/32 + 2/3/32 - 3/461 - 3/362 - 2/263, ... 

6„ = /3„ + ... + /3i/„ - 2/26„_i - ... - n/A, ... (5) 

so, by reparametrization, beyond the critical dimension (/3i 7^ 0) we can change any coefficient 
but /3i. We can fix any higher coefficient with zero value, if we take 

f _ n _ n2 n _ Pn + --- /^N 

J2 — -5-, /3 — 77-5- + /2 5 ... ; In — 7 TViT ^ ••• 

/3i 2/3i {n - l)/3i 

In this case we have exact classical dynamics in the (external) space-time and simple scale 
dynamics. 

We will consider in applications also the case when only one of the higher coefficients is 
nonzero. 

In the critical dimension of space-time /3i = and we can change by reparametrization 
any coefficient but /32 and P^. From the relations ([5]), we can define the minimal form of the 
RD equation 

A = P^A^ + (3^A\ (8) 

e.g. 64 = when 

/3 = ^ + ^/, + /|. (9) 

P2 P2 

/2 remains arbitrary and we can take e.g. /2 = 0. 
We can solve ([8]) as implicit function, 

uMh^-u ^ ^^ht^ ^ = 1 + (10) 

A P2 

than, as in the noncritical case, explicit solution will be given by reparametrization represen- 
tation ([3]). If we know somehow the coefficients /3„, e.g. for first several exact and for others 
asymptotic values (see e.g. [7]) than we can construct reparametrization function and find 
the dynamics of the running coupling constant. 



1.2 Renormdynamic functions (RDF) 

We call RDF functions Qn = fn{t), which are solutions of the RD motion equations 

9n = Pn{9), l<n<N. (11) 

In the simplest case of one coupling constant, the function g = f{t), is constant g = Qc when 
/3{gc) = 0, or is invertible (monotone). Indeed, 

9 = f{t) = r{f-\g))=m- (12) 

Each monotone interval ends by UV and IR fixed points and describes corresponding phase 
of the system. 

In the lattice (gauge) theory approach to the renormdynamics (see, e.g. [8]), recently 
running coupling constant dynamics were calculated for SU{2) Yang-Mills model [9]. The 
result is in agreement with perturbative calculations at small scales; at an intermediate scale 
the coupling constant reaches its maximum(~ 1.25); than decrease. So, at the maximum, we 
may have nontrivial zero of the /3— function, which corresponds to the conformal invariance 
of the gluodynamics at this point. Beyond this point we have another phase, strong coupling 
phase with decreasing coupling constant similar (identical?!) to the abelian (monopole?) 
theory. 

If we approximate the form of the curve near maximum as 

a{t) = ttc — b\t — tc 

for the /3— function we obtain 

a = /3{a,t) = sign{tc — t)bn{-^ — . (14) 

Of course this is not usual /5— function, function of a only. It depends also on t. For t > tcwe 
have perturbative phase. For n> 1, (3{ac,t) = 0. 

The fundamental quark and gluon degrees of freedom are the relevant ones at high temper- 
atures and/or densities. Since these degrees of freedom are confined in the low temperature 
and density regime there must be a quark and/or gluon (de) confinement phase transition. 

It is difficult to describe the phase transition because there is not known a local parameter 
which can be linked to confinement. We consider the fractal dimension of the hadronic/quark- 
gluon space as order parameter of (de) confinement phase transition. It has a value less than 
3 in the abelian, hadronic, phase, and more than 3, in nonabelian, quark-gluon, phase. 
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1.3 Hamiltonian extension of Renormdynamics 

The renormdynamic motion equations ( ITTi) can be presented as nonlinear part of a hamiltonian 
system with linear part 



= (15) 
09n 



the hamiltonian and canonical Poisson bracket are [TO] , 



N 



H=Y. /3(^)n^„, {9ni ^m} = ^nm (16) 



n=l 



In this extended version, we can define optimal control theory approach \TT\ to the unified 
field theories. We can start from the unified value of the coupling constant, e.g. a~^{M) = 
29.0... at the scale of unification M, in the minimal supersymmetric extension of the standard 
model (see e.g. [12]), put the aim to reach the SM scale with values of the coupling constants 
measured in experiments, and find optimal threshold corrections to the RD coefficients. 

Any Hamiltonian quantum (and classical) system can be described by infinitely divisible 
distributions because, in the functional integral formulation, we use the following step 

U{t) = e~''" = {e-'^"f (17) 

In the case of the scalar field theory, e.g. 

L(^) = Id.cpd^^ - - = g^^{^d,cj>d'^cj> - _ l^n)^ (18) 

so, to the constituent field (pN, when n > 2, corresponds higher value of the coupling constant, 

= gN"^ (19) 
For weak nonlinearity, n = 2 + 2e, gN = + eln + 0(e^)) 



2 Negative binomial distribution 

Negative binomial distribution (NBD) 

= 11^(1 - pYp" = EpL^i^^n^^r, i: = i. (20) 

1 [r)n\ 1 [r)nl < n > +r < n > +r 

provides a very good parametrization for multiplicity distributions in e~^e~ annihilation; in 
deep inelastic lepton scattering; in proton-proton collisions; in proton-nucleus scattering. 
Hadronic collisions at high energies (LHC) lead to charged multiplicity distributions whose 
shapes are well fitted by a single NBD in fixed intervals of central (pseudo)rapidity rj [13j . 
The generating function for NBD is 

m = (1 + ^^(1 - h)r = ( ^ " )^(l - Ph)r = E Pin)h- (21) 

An useful property of NBD with parameters < n > and r is that it is the distribution 
of a sum of r independent random variables with a Bose-Einstein distributioiiE] and mean 

< n > /r, 

< > +1 < n > +1 
^np{n) =< n >= _ , f{x) = ^xXn) = (1+ < n > (1 - x))"^ 

T = ^ , . , ~ /icj < n >, <n » 1. (22) 

Temperature defined in ( l22l) gives an estimation of the Glukvar temperature when it radiates 
hadrons. We see that universality of NBD in hadron-production is similar to the universality 



"'^A Bosc-Einstein, or geometrical, distribution is a thermal distribution for single state systems. 



of black body radiation. The generating function of a sum of independent random variables 
is the product of their generating functions. Indeed, for 

n = Til + n2 + ... + Hk, (23) 

with rii independent of each other, the probability distribution of n is 

P{n)= J2 5{n-Y.n,)p{ni)-p{nk). F{h) = Y.K''P{n) = f{hf (24) 

ni,...,nfc n 

An incoherent superposition of m emitters that have a negative binomial distribution with 
parameters r,< n > produces NBD with parameters mr,m < n >. So, 

F{r, < n >)™ = F{mr, m< n >) (25) 

We can put this equation in the closed nonlocal form 

^ ^ ™^ n ^ < n > d Xid Xod , . 

QgF = F\ Q, = q'', D = -+ =J- + ^ 26 

dr a < n > axi ax2 

Note that F{xi,X2) may be any function of the tj^e 

F{x,,x2) = fi^r = v{-r\ ^{x) = /(x)i (27) 

X2 X2 

By construction we know the solution of the nice equation (1261) as GF of NBD, F. We obtain 
corresponding differential equations, if we consider q = 1 + e, for small e, 

{D{D - 1)...{D - m + 1) - {lnFr)^> = 0, (-Ii^±l^ - (lnF)-)vl/ = 0, 

(D^ - $"^)vl/ = 0, m = l,2,3,..., D^ = ,^i^±lL,$ = inF, (28) 

L [U + 1 — m) 

with the solution = F = exp($). 

2.1 NBD, Poisson and Gauss distributions 

Fore high values of 0:2 = r the NBD distribution reduces to the Poisson distribution 

F(xi,X2,h) = (1 + —(1 - /i))-"2 e-"i(^-'*) = e-<">e''<"> = V P(n)/i", 
X2 

P(n) = e-<">— ^ (29) 
n! 



For the Poisson distribution 



\h=\ 



--< nin - 1) >=< n >' 



dK^ 

D"^ =<n^ > - <n >2=< n > . (30) 



In the case of NBD, we have the following dispersion low 



= - <n>^ + <n>, (31) 
r 



which coincides withe previous expression for high values of r. 



Poisson GF belongs to the class of the infinitely divisible distributions, 

F{h,<n>) = {F{h,<n> /k))'' (32) 
For high values of < n >, the Poisson distribution reduces to the Gauss distribution 

P n = e-<"> 1 — = , exp '-) 33 

Let us consider the values q = n,n = 1,2,3, ... and take sum of the corresponding equations 
we find 

^(-^)^ = r^-p (34) 

Now we invent a Hamiltonian H with spectrum corresponding to the set of nontrivial 
zeros of the zeta function, in correspondence with Riemann hypothesis, 

Tl Ti 

-Dn = 2 + ^ ^( 2 ^ ^" ^ ^^^^ ^ ^^^^ ^ "' ^ 

n 11 

= Hn= Yl Hi{xm), Hi{x) = i{- + xdx) = --{xp + px), p = -idx (35) 



m=l 



The Hamiltonian H = Hn is hermitian, its spectrum is real. The case n = 1 corresponds to 
the Riemann hypothesis. The case n = 2 corresponds to NBD, 

(^{l + iH^)F = -^, F(xi,X2;/i) = (1 + ^(1-/1))^^^ (36) 

I - F X2 

Let us scale X2 — )■ \x2 and take A — )■ oo in (l36l) . we obtain 

a\ + ^/7(x))e-(i-^)- = H{x) = ^(1 + x9.) = -^{xp + px), 

H{x)i)E = E^pE, i'E = cx~% s = ]^ + iE, c = 1/^271-, 

POO 

/ dxi)E{xyiJE>{x) = 5{E - E'), (37) 
Jo 

a-D)e-^ = al + ^H{x))e-'' = (38) 
dxx^-\{- + iH{x))e-^=<x^~'\a- + iH{x))e~^>= = r(s)C(s), 

< a;""'|C(^ + iH{x))e-'' >=< C(i - iH{x))x'-'\e-'' > 

= al - ^E) < x^-'\e-^ >= al - ^Ens), C(^ - ^E) = C{s). (39) 



A slightly different consideration is the following. If we rescale x ^ xy in (l38ll . multiply 
by y^~^ and integrate by y, we obtain usual integral formula for zeta-function 



s-l 



a-D) / y^-'e-^ydy = / dy^-- 
Jo Jo e^y — 1 

a-D)x-T{s) = x-^ dy^^--, 

Jo ey — 1 

1 Is) JO ey — 1 



3 Renormdynamical formulation of z - Scaling 



In the z - Scaling approach to the inclusive multiparticle distributions (MPD) (see, e.g. [H]), 
different inclusive distributions depending on the variables described by universal 

function '^{z) of a fractal variable z = It is interesting to find a dynamical system 

which generates these distributions and describes corresponding MPD. 

We can find a good function if we know its derivative. Let us consider the following RD 
like equation 

^ = V{%, / -— = In -, V{^) = -P^{z) + 7*(^)'+^ (41) 

CtZ "'*(^o) V[Xj Zq 

Corresponding solution for \1/ is 

^(z) = xl>(z;/3,7,c,n) = (1 + cz^f^)-^ (42) 

where the integration parameter c is defined from the normalization condition on \1/ as 

c^- = (^)'^ ^ (43) 
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Figure 1: Typical form of z-scaling distribution ( l42l) . ^{z, 9, 9, 1, 1) 



4 Multiparticle production stochastic dynamics 

Let us imagine space-time development of the the multiparticle process and try to describe 
it by some (phenomenological) dynamical equation. We start to find the equation for the 
Poisson distribution and than naturally extend them for the NBD case. 

Let us define an integer valued variable n{t) as a number of events (produced particles) 
at the time t, n(0) = 0. The probability of event n(t), P{t,n), is defined from the following 
motion equation 

P, ^ = r(P(t, n-1)- P{t, n)),n>l 



Pi(t,0)) = -rP(t,0), 
P(t,n) = 0, n < 0, 



(44) 



so 

P{t,0) = Po{t) = e--\ 
P{t,n) = Q{t,n)Po{t), 

Qt{t, n) = rQ{t, n-1), Q{t, 0) = 1. (45) 
To solve the equation for Q, we invent its generating function 

F(t,/i) = ^/i"Q(t,n), (46) 

n>0 

and solve corresponding equation 

F, = rhF, F{t, h) = e^'^ = ^ /i"^, Q(t, n) = (47) 



so 



-Art) 



P(t,n) = e-^*^ (48) 

is the Poisson distribution. If we compare this distribution with ( l33l) . we identify < n >= rt, 
as if we have a free particle motion with velocity r and the distance is the mean multiplicity. 
This way we have a connection between n-dimension of the multiplicity [15] and the usual 
dimension of trajectory. As the equation gives right solution, its generalization may give more 
general distribution, so we will generalize the equation fH4|) . For this, we put the equation in 
the closed form 

Pt{t,n) = r(e-^" - l)P{t,n) 

= J2D,d'P{t,n), D, = {-lf^, (49) 
fc>i 

where the D^, k > 1, are generalized diffusion coefficients. 

For other values of the coefficients, we will have other distributions. 

4.1 Fractal dimension of the multiparticle production trajectories 

For mean square deviation of the trajectory we have 

<{x- >=< > - <x >'^= D{xf ~ t^/d/^ (50) 

where df is fractal dimension. For smooth classical trajectory of particles we have dj = 1] 
for free stochastic, Brownian, trajectory, all diffusion coefficients are zero but D21 we have 
df = 2.\a the case of Poisson process we have, 

D{nf =<n^ > ~ <n t, dj = 2. (51) 

In the case of the NBD and KNO distributions 

D{nf ~ t^, df = 1. (52) 



As we have seen, rasing k, KNO reduce to the Poisson, so we have a dimensional (phase) 
transition from the phase with dimension 1 to the phase with dimension 2. It is interesting, if 
somehow this phase transition is connected to the other phase transitions in strong interaction 
processes. 

For the Poisson distribution GF is solution of the following equation, 

F = -r(l - h)F, (53) 
For the NBD corresponding equation is 

If we change the time variable as t = T'^-f , we reduce the dispersion low from general fractal 
to the NBD hke case. Corresponding transformation for the evolution equation is 

Ft = -dfT'^f-^R{T'^^)F, (55) 

we ask that this equation coincides with NBD motion equation, and define rate function R{T) 

dfT^^f-'RiT'^n = (56) 

now the following equation defines a production processes with fractal dimension dp 

F, = -R{t)F, R{t) = ~ (57) 

It is a pleasure to thank Yu.M.Bystritskiy, M.V.Tokarev and the members of the seminar 
QNMtt for stimulating discussions and various help. 
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